Fluctuations assisted stationary entanglement in driven quantum systems 
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We analyze the possible quantum correlations between two coupled dimer systems in the presence 
of independent losses and driven by a fluctuating field. For the case of the interaction being of a 
Heisenberg exchange type, we first analytically show the possibility for robust stationary entangle- 
ment for realistic values of the dissipation rates and then analyze its robustness as a function of the 
noise to signal ratio of the pump. We find that for a common fluctuating driving field, a stochastic 
resonance effect appears as function of the ratio between field strength and noise strength. The 
effect disappears in the case of uncorrelated or separate pumps. Our result is general and could 
be applied to different quantum systems ranging from electron spins in solid state, to ions trap 
technologies and cold atom set ups. 
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Introduction: Entanglement is nowadays commonly 
considered a resource for quantum information process- 
ing [l] . Being of a purely quantum nature can be easily 
degraded by reservoir contaminations. Nevertheless, re- 
cently the possibility of having stationary entanglement 
in open quantum systems has been put forward p]-[5] , also 
achieving experimental evidence [6 . However it is not 
clear how robust is such stationary entanglement against 
other sources of noise that can eventually affect the sys- 
tem, especially when originating from the pump. 

In this work we aim at clarifying the robustness 
of stationary entanglement in lossy coupled systems 
in the presence of noise in the amplitude fluctuations 
of the driving field. We assume the case of the al- 
ways on exchange interactions between the systems, 
usually employed when dealing with local dissipative 
environments [4 . The entanglement robustness is char- 
acterized by means of a figure of merit related to the 
quantum concurrence and accounts for the signal to noise 
ratio of the driving field. We distinguish the case of lo- 
cal and global fluctuations and show that in the latter 
case stochastic resonance like effects [^ can emerge for 
an optimal ratio between the strength and the noise of 
the common driving field. The effect disappears in the 
case of uncorrelated or separate pumps. 

Our analysis is general, the results are analytic, and 
could be applied to different set ups ranging from electron 
spins in double quantum dots or in effective spin models 
generated in driven ions, cold atoms and coupled cavity 
arrays upsfSUTT] 

The system: We assume two coupled dimer systems 
labelled by A and B interacting through an exchange in- 
teraction. The systems are subject to incoherent driving 



of amplitude a from outside and dissipate independently 
into two separate reservoirs with equal dissipation rates 
given by 7 (see Fig. 1). The Hamiltonian describing the 
coherent part of the interaction in this case is: 



H = J{(j\(JB + CFACT^b) + ^(^A + ^b)' 



(1) 



with as the usual Pauli operators on C^ such that a = 
(cr^ + i(7y)/2 and 2crtcr - / = c^^ 

The open system dynamics will be described by the 
master equation 



where 



i[H,p\^V[cja]p^V[cjb]p, (2) 



V[a\h := '^{2aho) — o) ah — ho/ a)^ (3) 



denotes the dissipative super-operator and 7 is the decay 
rate into the two separate environments. 

In this work besides looking for entanglement at steady 
state, we would also investigate its robustness against 
fluctuations in the driving field. We thus assume that the 
driving field on top of its coherent amplitude a, exhibits 
a Gaussian white noise term ^(t), with 



and ?7 > measuring the noise strength. 
Hence the master equation Q will become 



(4) 
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p = -i[H, p] + V[aA]p + V[(Jb]p + M[(t\- (jI]p, (5) 

with N the superoperator describing the noisy effects of 

The steady state solution of Eq.([5| can be 
found by writing the density operator and the 
other operators in a matrix form, in the basis 



I^Common noisy driving 



Independent noisy driving^ 




FIG. 1: The system under consideration. Two coupled dimer 
systems (atoms or quantum dots or spins) are interacting 
and dissipating in different reservoirs. In b) two uncorrelated 
noisy pumps are continuously driving the systems whereas in 
a) the same pump is driving both. 



B = {|0)a|0)b,|0)a|1)b,|1)a|0)b,|1)a|1)b}. Let us 
parametrize the density operator as 



a 61 + ib2 c\ + IC2 d\ 

bi - ib2 e /i + 1/2 gi 

ci - iC2 /i - if2 h ii 

di — id2 gi — ig2 ii — ii2 1 — « 




(6) 



Once we know the steady state density matrix, to 
quantify the entanglement we will use the concurrence 



m 



C{pss) ■= max {0, Ai - A2 - A3 - A4} , 



(7) 



where A^'s are, in decreasing order, the nonnegative 
square roots of the moduli of the eigenvalues of PssPss 
with 



K®4)pLK®4). 



(8) 



and p*g denotes the complex conjugate of pss- 

We will separate our analysis in two cases, the common 
driving case where both systems are subject to the same 
fluctuating pump and the case when independent driving 
is applied to each one. 

Common driving: If the driving field is common to 
both systems (as also in |13j), then the two systems ex- 
perience the same fluctuations. Hence the superoperator 
A/" reads (see also [14]) 

■A/"!^^; cr^jp = -V [K + ^1). [K + <l p]] ■ (9) 

We can then analytically solve the linear system of equa- 
tions coming from the master equation ([5| 

= -iJ[{a\aB +crAO-J^),p] - ia[{a\^al),p] 
-V [K + ^1). [{^'a + ^1). P]] + ^[^a]p + V[ctb]p. (10) 

for a,6i,62,ci,C2,c^i,6^2,e,/i,/2,^i,^2,^,n,^2. The ex- 
plicit solution is reported in Appendix A. In Figure |2] we 




FIG. 2: Steady state entanglement as quantified by concur- 
rence C(a, 77) vs a and 77 for J = 2 in units of the dissipation 
rate 7. 




FIG. 3: The quantity A (a, 77) vs a and 77 for J — 2 (units as 
in Fie;. 2). 



plot the concurrence C [12 vs amplitude of the pump a 
and fluctuation 77 for a fixed nonzero values of J. 

We note here that C is not always monotonically de- 
creasing vs 77 which can been seen by introducing the 
quantity 



A(a,7?)=C(a,r,)-C(a,7? = 0) 



(11) 



i 



for C(a, 77) — C(a, 77 = 0) > 0, which quantifies the 
difference (if positive) between the concurrence in pres- 
ence and absence of fluctuations in the driving field. 
A(a,7/) = when C{a,r]) - C{a,r] = 0) < The in- 
verse U shape of such a quantity vs 7^, as shown in Fig. 
^ indicates a stochastic resonance-like effect [7 . That 
is, there is an optimal non zero value of noise strength 
for the driving field which maximize the concurrence. In 
Fig. [2] we may also observe that entanglement is quite 
robust in the range of pumping slightly larger than the 
decay rate and it survives for white noise up to 10% in 
the fluctuating field. 

To quantify the robustness of entanglement versus 
noise in the driving field we may consider the signal 
to noise ratio SNR := a/y^ for a corresponding to 
maXcK C{a^r] = 0) and rj corresponding to the value for 
which maxo, C(a, 77) becomes zero. The smaller is this 
quantity, the more robust is entanglement. The SNR 
has been plotted in Figure El (bottom curve) vs J. The 
value J = 1 is optimal in the sense that admits a larger 




FIG. 4: The quantity SNR is plotted vs J in units of the 
dissipation rate 7. Top (bottom) curve refers to the separable 
(non separable) case. 




FIG. 5: Concurrence C{a,ri) vs a and 77 for J — 2. 



interval of r] values for which maxQ, C{a^r]) > 0. 

Independent driving fields: If the driving field is sepa- 
rate for the two systems, then they will experience inde- 
pendent fluctuations. Hence the superoperator A/" reads 
(see also [14 ) 

■A/"!^^; ^|]p = -V [<^\, Wa. pW -vK, [^I, p]] ■ (12) 

We can then again analytically solve the linear system of 
equations coming from the corresponding master equa- 
tion for the steady state 

= -iJ[{a\aB ^ CTAcr]^), p] - la [{a\ ^ a%) , p] 
-V W% W% P]] -VK. K.P]] + V[aA]p + ViasMU) 

for a,bi,b2,Ci,C2,di,d2,eJij2,gi,g2,h,ii,i2 (the solu- 
tion is reported in Appendix B) and fully characterize 



the state of the system. 

In Figure |5] we plot as before the concurrence vs the 
pump strenthg a and noise 77 for a fixed nonzero values of 
J. In this case C is a monotonically decreasing function 
of T]. Hence, we do not have any stochastic resonance 
like effects, as matter of fact the quantity (11) is zero 
everywhere for all values of signal to noise ratio. 

By comparing Figure [5] and |2] we may realize that en- 
tanglement is less robust in this case with respect to the 
previous one. As matter of fact, the maximum value of 
concurrence decreases to zero when the signal-to-noise 
ratio is ^ 1.8 in contrast to the previous case where the 
corresponding ratio was ^ 1.4. In Figure [I] (top curve) 
we have plotted the quantity SNR vs J. Likewise the 
previous case the optimal value is J = 1. 

Our results in general are in agreement with similar 
phenomena that arise in spin chains from the interplay 
of dissipative and dephasing noise sources [15 and to 
those introduced in Ref.[3]. In the latter, two dissipat- 
ing cavity modes interacting with an atom were driven 
by thermal noise. However, also in this case the additive 
noise besides the local dissipative fluctuations is common 
to both systems, hence we can quite generally argue that 
stochastic resonance effects on entanglement arise from 
the interplay of separate and common fluctuations of dif- 
ferent kind. Moreover, we can speculate that non-white 
noise in Eq.Q could lead to more pronounced stochas- 
tic resonance like effects and would like to study this in 
follow ups works. 

Conclusion: We have studied the robustness of entan- 
glement that can be created in the steady state of a driven 
two coupled dimers interacting with an exchange inter- 
action. As a further noise source we have considered 
amplitude fluctuations of the driving field. The entan- 
glement robustness has been characterized by means of a 
figure of merit accounting for the signal to noise ratio of 
the driving field. We have distinguished the case of local 
and global fluctuations and showed that in the latter case 
stochastic resonance like effects can emerges. Our anal- 
ysis could be applied to effective exchange interaction 
models between dimers realizable in various implemen- 
tations in solid state, ions traps, cold atoms and Cavity 
QED set ups. 



Appendix A: Solution of Eq.(lO) 



D = 24a^ + J^{1 + 2r]f{l + 6r]) + 4a^ (l5 + 82r] + lUrj'^) + 2{a + 2^7^)^ (21 + 184r] + 4327^^) 
+(1 + 2r]f (9 + 1507? + 8807?^ + 20807?^ + 15367?"^) 



+2J2(1 + 67?) {2a^ + (1 + 27?)^ (5 + I67? + 87?^) + a^ (5 + I87? + 167?^)) 



Da 



6a^ 



a^{9^ 66r] + 1447?^ + J^(l + 67?)) + 37?^ (1 + 27?) 



Dbi 

Db2 

ci 
Ddi 

Dd2 

De 



X (9 + J"^ + 967? + 3047?^ + 2567?^ + 2J^ (5 + I67? + 87?^)) 
+2a^7? (J^ (7 + 237? + 247?^) + 3 (9 + 637? + 1527?^ + 1447?^)) 
= a {Ga"^ + a^ (9 + 487? + 1207?^ + J^(l + 67?)) + 7?(1 + 27?) (J^(ll + 127?) + 3 (9 + 6O7? + 647?^))) 



jQf7? (3 

bi 



h 

Dgi 

Dg2 

h 

H 
i2 



lOa^ + 227? + 327?^ + J^(3 + 67?)) 



= a^(io + 87?) + a^ (9 + 487? + 1207?^ + 1607?^ + J^ (l + 47? + 87?^)) 
+7?(1 + 27?) (9 + J"^ + 967? + 3047?^ + 2567?^ + 2J^ (5 + I67? + 87?^)) 

= -Ja^ (9 + 6a^ + 667? + 967?^ + J^(l + 67?)) 

= 6a^ + a^ (15 + 747? + 1447?^ + J^(l + 67?)) + 7? (l + 57? + 67?^) 
X (9 + J"^ + 967? + 3047?^ + 2567?^ + 27^ (5 + I67? + 87?^)) 



W (9 + 1027? + 4987?^ + 10727?^ + 8647?"^ + J^ (l + I87? + 547?^ + 487?^)) 

di 



a (9 + Ga^ + 1057? + 4187?^ + 6807?^ + 3847?"^ + 5a^ (3 - 

+aJ^ (1 + 137? + 347?^ + 247?^ + a^(l + 67?)) 



h 167? + 247?^)) 
-Ja (9 + 8I7? + 2067?^ + 1607?^ + a^(6 + 227?) + J^ (l + 57? + 67?^)) 



e 

92 



(Al) 



(A2) 
(A3) 
(A4) 
(A5) 
(A6) 

(A7) 
(A8) 



(A9) 
(AlO) 
(All) 

(A12) 
(A13) 
(A14) 
(A15) 
(A16) 



Appendix B: Solution of Eq.(13) 



D = J^(l + 27^)^(1 + 4r]) + (3 + 8r]) (l + 2a^ + 6/? + 8r]^)^ (3 + 2a^ + lOr] + 8r]^ 



Da 



Dbi = 

Db2 = 

ci = 

C2 = 

Ddi = 

Dd2 = 

De = 



h 

Dgi 

Dg2 

h 

h 
i2 



+2J2(1 + 2r]) {a^{2 + 8r]) + (1 + 2r]f (5 + 32r] + 56/?^ + 32r]^) + a^ (5 + 38/? + 96r]^ + 64/^^)) (Bl) 

2a^(3 + 8r]) + a"^ (9 + 66/? + 184/^^ + 192/^^ + J^ (l + 6/? + 87?^)) 

+7?^ (1 + 67? + 87?^) (9 + J"^ + 727? + 1767?^ + 1287?^ + 27^ (5 + 127? + 87?^)) 

+2a^7? (9 + 937? + 3527?^ + 5927?^ + 3847?"^ + J^ (3 + 237? + 487?^ + 327?^)) (B2) 

a (2a'^(3 + 87?) + 7? (l + 67? + 87?^) (9 + 367? + 327?^ + 7^(5 + 47?))) 

+a^ (9 + 6O7? + 1447?^ + 1287?^ + J^ (l + 67? + 87?^)) (B3) 

Ja7? (J^ (1 + 67? + 87?^) - (3 + 87?) (1 + 2a^ + 67? + 87?^)) (B4) 

bl (B5) 

62 (B6) 

a^ {J\l + 27?) + (3 + 87?) (3 + 2a^ + IO7? + 87?^)) (B7) 

- {Ja^ (J' (1 + 67? + 87?') + (3 + 87?) (3 + 2a' + 147? + 87?'))) (B8) 
2a^(3 + 87?) + ^"^(1 + 27?) (15 + 767? + 967?' + j'(l + 47?)) 

+7? (1 + 77? + 147?' + 87?^) (9 + J"^ + 727? + 1767?' + 1287?^ + 2j' (5 + 127? + 87?')) 

+a' (9 + 1147? + 5707?' + 14087?^ + 16967?"^ + 7687?^ + J' (l + I87? + 787?' + 1287?^ + 647?"^)) (B9) 

di (BIO) 

(Bll) 

a (1 + Of' + 57? + 47?') (j' (1 + 67? + 87?') + (3 + 87?) (3 + 2a' + IO7? + 87?')) (B12) 

- {Ja (j' (1 + 77? + 147?' + 87?^) + (3 + 87?) (a' (2 + 67?) + 3 (l + 77? + 147?' + 87?^)))) (B13) 
e (B14) 

91 (B15) 

92 (B16) 
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